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Abstract 

We explore a method developed in statistical physics which has been argued to have expo¬ 
nentially small finite-volume effects, in order to determine the critical temperature T c of pure 
SU(3) gauge theory close to the continuum limit. The method allows us to estimate the 
critical coupling /3 C of the Wilson action for temporal extents up to N r ~ 20 with <0.1% 
uncertainties. Making use of the scale setting parameters r 0 and in the same range of ,0- 
values, these results lead to the independent continuum extrapolations T c r 0 = 0.7457(45) and 
Tcy/to = 0.2489(14), with the latter originating from a more convincing fit. Inserting a conver¬ 
sion of ro from literature (unfortunately with much larger errors) yields T c /Am§- = 1.24(10). 
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1. Introduction 


Even though light quarks play an essential role for the phenomenological understanding of 
heavy ion collision experiments it can be argued that, due to their large multiplicity in the 
initial state and their Bose-enhanced distribution functions in the plasma phase, gluons are 
the single most important degree of freedom influencing the formation and evolution of QCD 
matter. Gluons are also much easier to study with non-perturbative lattice methods than 
light quarks. Therefore studies of pure SU(3) gauge theory at high temperature continue 
to constitute an important laboratory system, both for developing numerical techniques and 
for gaining physics understanding on observables where a high precision is needed. Recent 
examples of topics studied include scale setting, renormalization, and methods for statistical 
error reduction (cf. e.g. refs. [1]- [5]). Our own interest stems from attempts to measure 
real-time observables such as transport coefficients [6]- [8], in which case theoretically well- 
founded methods [9] can probably be applied (if at all) only after the infinite volume and 
continuum limits have been reached with a high precision [10]. 

In the present contribution we use the pure SU(3) gauge theory as a test bench for studying 
finite-volume scaling in the vicinity of a first-order phase transition. Concretely, our primary 
goal is to determine the critical coupling /3 C for values of N T much larger than have been 
achieved before (here N T = l/(aT) is the number of lattice points in the periodic imaginary- 
time direction; a is the lattice spacing; and T is the temperature). Let us remark that values 
of (3 C as a function of N r have attracted recent interest as tests of semi-analytic models [11,12], 
and indeed new high-precision values at large N r put the functional dependences predicted 
by these frameworks under tension [7]. 

The second focus point of our study is that of scale setting [13]. In particular, we consider 
two scales that have been frequently employed, denoted by r 0 [14] and [15]. Neither of 
these scales has a direct physics interpretation; however they are relatively straightforward 
to measure and can in principle be related to physical quantities in a separate study. On 
the other hand, in the thermal context there is one directly physical quantity, the critical 
temperature T c , which would have certain advantages as a scale setting parameter, permitting 
for instance for an easy comparison of theories with different matter contents but with similar 
macroscopic properties (this assumes, of course, that all theories considered have a sharply 
defined transition point). Therefore, we make use of our results in order to obtain a largely 
independent estimate for T c r 0 [16] and a new estimate for T cy /t^. It should be acknowledged, 
however, that close to the continuum limit we also see indications of growing systematic 
uncertainties, particularly in case of r 0 . 

The plan of this note is the following. After introducing and testing the basic method of 
our study in sec. 2, we employ it in order to estimate the critical coupling /3 C as a function of 
N t in sec. 3. The issue of scale setting is addressed in sec. 4, and we conclude in sec. 5. 
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Figure 1: Left: Determination of the right-most minimum (open circle) from the distribution 
of ReP. Middle: The corresponding triangle separating the two phases, with the vertical 
line placed at the position of the open circle. Right: The resulting function s(/3) (cf. eq. (3)), 
permitting for an estimate of /3 C from the crossing of zero. The statistics of each data point 
is 0(1O 5 ) sweeps; statistical errors are based on jackknife estimates. 


2. Method 


The Wilson plaquette action, 

5 w^E Tr (l- P /-)’ (!) 

studied on an N r x Nf lattice with periodic boundary conditions in all directions, has a global 
Z(3) symmetry that is broken at and above the transition point for N s —> oo. We denote the 
location of the transition point by f3 c . Theoretical arguments [17] and empirical evidence [18] 
suggest that this is a first order phase transition. 

It has been shown through a study of g-state Potts models in three dimensions [19,20] that 
even though most observables, such as susceptibilities, show powerlike finite-volume effects 
at a first-order transition point, there is a particular definition of a pseudocritical point for 
which finite-volume effects are exponentially suppressed. This is obtained if the “weights” of 
the phases with no degeneracy ( w c ) and with g-fold degeneracy (w d ) are related through 

qw c = w d . (2) 


The weight can be defined through the “volume” of the distribution of some observable which 
has a good overlap with the order parameter. More formally, the weight corresponds to the 
partition function associated with the phase considered. 

For SU(3), a suitable observable is the Polyakov loop expectation value. Carrying out 
measurements in the vicinity of /3 C , we define 


3 w c - w d 

s(P) = x--- ■ 

3 w c + w d 


(3) 
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Figure 2: The pseudocritical couplings extracted from our method at N r = 4 (closed circles), 
normalized to the central value of the infinite-volume estimate /3 C = 5.69254(24) from ref. [18]. 
We also compare with susceptibility maxima from ref. [21] (open squares). The grey band 
illustrates our infinite-volume extrapolation (constant fit to N S /N T > 3). 

By construction s(/3) equals +1 deep in the confined and —1 deep in the deconfined phase. 
The critical point is obtained by interpolating to the location where s(/3 c ) = 0. 

In order to implement the idea, we need to introduce a criterion for separating a distribution 
into contributions from different phases. In a finite volume, when the distributions overlap, 
the procedure is not unique. In this study, we define a separatrix by looking for a minimum in 
the distribution of Re P, where P denotes the Polyakov loop (cf. fig. 1 (left)). This minimum is 
employed for defining a triangle separating the two phases (cf. fig. 1 (middle)). The resulting 
weights are the inputs for eq. (3); /3 C is obtained by a linear interpolation from points on both 
sides of the zero (cf. fig. 1 (right)). 

The results obtained with this procedure are shown in fig. 2 for N T = 4. They have been 
normalized to a classic value from ref. [18], and are compared with recent high-precision 
pseudocritical points extracted from Polyakov loop susceptibility maxima [21]. We conclude 
that for N s > 3 N r , no finite-volume effects can be observed within our resolution (~ 0.005%). 
For N s < 3N t1 we expect j3 c to be slightly underestimated. 
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3. Results at finite lattice spacing 

We carried out measurements for 4 < N r < 22, increasing N r in steps of two. We computed 
on several volumes for ensembles with N T < 18, verifying that volume dependence is below 
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Nr 

Pc [18,22] 

Pc [21] 

Pc [23] 

f3 c [our value] 

N s used 

-^total 

4 

5.69254(24) 

5.692469(42) 


5.69275(28) 

14,...,40 

83xl0 6 

5 



5.8000(5) 




6 

5.8941(5) 

5.89410(11) 


5.89425(29) 

20,...,40 

28x 10 6 

8 

6.0624(10) 

6.06212(44) 


6.06239(38) 

28,32 

4.2x 10 6 

10 




6.20873(47) 

32,...,56 

15x 10 6 

12 

6.3380(17) 



6.33514(45) 

CM 

F- 

o' 

21x 10 6 

14 




6.4473(18) 

48,56 

12x 10 6 

16 




6.5457(40) 

64 

2.5xl0 6 

18 




6.6331(20) 

56,64 

3.6x10® 

20 




6.7132(26) 

64 

4.0x10® 

22 




6.7986(65) 

64 

5.9x10® 


Table 1: The infinite-volume critical points of SU(3) gauge theory according to various stud¬ 
ies. lV tota i indicates the total numbers of configurations (all volumes and values of f3). Our 
data are based on constant fits to N s > 3 N r whenever several volumes are available. For 
N r = 4, 6 we have artificially enlarged the errors to account for systematics related to expo¬ 
nentially small volume corrections (cf. the text). 


statistical uncertainties. Subsequently we fit the data at N s > 3 N r to a constant. Given the 
resources at our disposal we used a single spatial extent N s = 64 for N T = 20, 22. Here, minor 
finite-volume effects start to contaminate our results. However, based on fig. 2, we expect 
the effects from a simulation with N s /N r = 64/22 = 2.9 to be below the 0.01% level, thereby 
being much below statistical errors. In contrast, at the smallest N r where statistical errors 
are extremely small, we have artificially saturated the errors at a constant value ~ 0.005%, 
corresponding to the expected uncertainty from finite-volume effects. Our final results at 
fixed N t , together with previous estimates from the literature, are collected in table 1. 


4. Continuum extrapolations 

In this section we convert the lattice-specific numbers of table 1 to values of T c in physical 
units. In order to achieve this two different scale setting parameters are considered, r 0 and 
y/t//, with the latter leading to a noticeably better description of the thermal data (cf. sec. 4.2). 

4.1. Scale r 0 

The scale ro/a [14] has been measured as a function of (3 in refs. [24,25] (see ref. [26] and ref¬ 
erences therein for previous work). We complement these results by a new set of simulations, 
with parameter values and results listed in table 2. The measurements were separated by 
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p 

r 0 /a [24] 

r 0 /a [25] 

r 0 /a [our value] 

N t x IV 3 

Aconf 

5.7 

2.922(9) 





5.8 

3.673(5) 





5.95 

4.898(12) 





6.07 

6.033(17) 





6.2 

7.380(26) 





6.3 



8.52(4) 

32 x 32 3 

216 

6.3 



8.51(2) 

32 x 48 3 

211 

6.3 



8.52(2)* 

32 x 64 3 

202 

6.336 



8.95(3) 

64 x 32 3 

220 

6.4 



9.80(3) 

36 x 36 3 

206 

6.5 



11.16(2) 

44 x 44 3 

202 

6.57 


12.18(10)** 




6.69 


14.20(12)** 




6.81 


16.54(12)** 




6.92 


19.13(15)** 





Table 2: The results for r 0 /a that have been used in our analysis. For (3 = 6.3 only the 
largest volume (indicated with an asterisk) has been included in subsequent fits. The values 
from ref. [25], marked with a double asterisk, do not come directly from r 0 but rather another 
scale r c , which has been converted into r Q through a continuum relation, whose systematic 
uncertainties are included in the errors. 


500 heatbath-overrelaxation updates. A number of standard techniques for statistical error 
reduction [27-29] were implemented in order to obtain these results. The static potential 
is extracted from Wilson loops with an ansatz based on two exponentials. The distance 
appearing in the static potential is tree-level improved [25], and subsequently a B-spline in¬ 
terpolation is carried out in order to extract r 0 /a from its definition [14]. (Note that due to 
the several steps involved, measurements are costly and systematic errors are difficult to get 
fully under control, particularly at large /?.) 

In order to permit for a subsequent interpolation, our data and older values [24,25] are fit 
in the range P € (5.7,6.92) to a rational ansatz inspired by ref. [30]: 


In 



P b\ . /66 0 \1 1 + ci/P + c 2 /P 2 

126 0 2 bl 11 V P J 1 + c 3 /P + c 4 /P 2 ’ 


(4) 


where bg = ll/(47r) 2 and b\ = 102/(47r) 4 . The fit parameters obtained read 1 


d = -8.17273 , c 2 = 14.9600 , c 3 = -3.95983 , c 4 = -5.30334 , x 2 /d.o.f. = 0.7 . 

(5) 

1 For the sake of reproducibility of subsequent results we show more digits than are statistically significant. 
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p 

(t 0 /a 2 ) wilson 

(f 0 /a 2 ) Wilson imp ' 

(f 0 /a 2 ) clover 

N t x IV 3 

A^conf 

5.6923 

0.6109(10) 

0.8234(9) 

1.0124(11) 

16 x 16 3 

455 

5.6923 

0.6103(7) 

0.8229(6) 

1.0119(7) 

16 x 24 3 

313 

5.6923 

0.6095(5) 

0.8220(5) 

1.0104(6) 

16 x 32 3 

248 

5.6923 

0.6010(4) 

0.8226(4) 

0.9905(4) 

24 x 32 3 

233 

5.6923* 

0.6097(3) 

0.8223(3) 

0.9800(4) 

32 x 32 3 

221 

5.8941 

1.9520(22) 

2.0989(22) 

2.2889(24) 

24 x 24 3 

465 

6.0625 

3.7129(39) 

3.8507(39) 

4.0626(41) 

32 x 32 3 

673 

6.2083 

5.9521(65) 

6.0873(66) 

6.3284(68) 

40 x 40 3 

476 

6.3352 

8.668(11) 

8.802(11) 

9.076(12) 

48 x 48 3 

315 

6.4487 

11.958(18) 

12.091(18) 

12.397(18) 

56 x 56 3 

254 

6.5509 

15.769(23) 

15.901(23) 

16.240(24) 

64 x 64 3 

305 

6.7130 

24.222(35) 

24.355(35) 

24.752(36) 

80 x 80 3 

250 


Table 3: Our results for t 0 /a 2 . The /3-values correspond approximately to those in table 1 
(apart from N r = 18, 22), with N r scaled up by a factor 4 in each case. For (3 = 5.6923 only 
the largest volume (indicated with an asterisk) has been included in subsequent fits. 


Based on the above equation, we convert the results in table 1 to values of r 0 T c : r 0 T c = 
(r 0 /a){f3 c )/N T . Subsequently we perform the extrapolation (a/r 0 ) 2 —l 0 using a fit quadratic 
in (a/r 0 ) 2 , illustrated in fig. 3(left), with the result 

r 0 T c = 0.7457(45) , x 2 /d.o.f. = 6.7 . (6) 

The error includes a rough estimate of systematic effects, encompassing the central values 
obtained by replacing the representation in eq. (4) through ln(r 0 /a) = X}n=o a n{P~ 6.25) n ; by 
carrying out the continuum extrapolation with a cubic fit; and by omitting f3 c corresponding 
to N t = 4. The first method increases the central value (T c r 0 ~ 0.7496), the second and 
third decrease it (T c r 0 ~ 0.7412,0.7424, respectively). However, in the first case the quality 
of the continuum fit decreases further from the already poor one in eq. (6), whereas in the 
second case the scatter of the data in fig. 3(left) suggests that including too much freedom 
in the fit distorts the outcome. A possible reason for the poor description of the data close 
to the continuum limit could be that estimates of r Q /a at (3 > 6.4 are systematically on the 
low side (by ~ 0(1%)), but unfortunately we have not been able to confirm this suspicion. 

The result in eq. (6) can be compared with r 0 T c ~ 0.7470(7) obtained in ref. [7] based on 
peak positions of Polyakov loop susceptibilities (here only statistical errors were included) 2 , 
as well as with the earlier value r 0 T c = 0.7498(50) [16]. 

2 For fixed N T the results of ref. [7] are consistent with the present ones, however their uncertainties from 
finite-volume effects are larger and only values up to N T = 16 could be reached. Therefore systematic errors 
would be larger than in the present study (but are more difficult to estimate reliably). 
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Figure 3: Left: Continuum extrapolation of T c r 0 , based on the data in table 1 and the 
interpolation from eq. (5). Right: Analogous results for T cy /t^ from the “Wilson” (open 
circles) and “Wilson imp.” (closed circles) discretizations, interpolated according to eq. (7). 

Finally, we recall that e.g. the values r 0 Aj^g = 0.586(48) [25], t-qA^s = 0.602(48) [31], 
r o^-Ms = 0.614(6) [32], and r 0 A^g = 0.637(32) [33] can be found in the literature (the third 
relies on the applicability of tadpole-improved lattice perturbation theory and the fourth of 
continuum perturbation theory at hadronic scales). Using the second value yields T c /A^g = 
1.24(10). Unfortunately the error is dominated by that in the relation of r 0 and A^g, so our 
new result in eq. (6) does not help to improve on previous estimates. 

4.2. Scale yTjj 

The scale is defined through the time that it takes for Wilson flow to adjust a chosen 
observable (= E ) to a predefined value [15]. We measured t 0 for a number of /? ~ j3 c , as listed 
in table 1. To study possible systematic effects, we made use of three different implementa¬ 
tions of E, based on the standard plaquette, tree-level improved, and clover discretizations, 
all of which are available within the DD-HMC package [34], Like for r 0 , the measurements 
were separated by 500 heatbath-overrelaxation updates; the volumes and the numbers of 
configurations used for measurements are shown in table 3. 

Given that the (3 values of table 1 correspond to the critical point, a set of fixed physical 
volumes can be chosen by scaling the corresponding N r by a constant amount. Setting 
N s = 4 N t we ensure that the box size is L = 4/T c ~ 5.3r 0 . For the smallest (3 we have 
carried out test simulations also at larger volumes, finding consistent results apart from the 
“clover” discretization for which volume dependence on the 3% level is visible. For our final 
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results we quote only those obtained with the two variants of the “Wilson” discretization that 
did not exhibit any volume dependence within statistical precision. Nevertheless systematic 
errors do grow with /?, because a longer integration trajectory in t is needed and because 
autocorrelation times tend to grow. 

As before, we represent the data as in eq. (4) for the interpolation, only this time replacing 
?’o —> y/t^- The resulting parameters are (for the “Wilson imp.” discretization) 3 

ci = -10.2116 , c 2 = 25.6819 , c 3 = -5.64462 , c 4 = 2.26845 , xVd.o.f. = 2.3 . (7) 

With this interpolation the critical values in table 1 can be converted into T CA /t^\ results are 
shown in fig. 3(right). A fit quadratic in a 2 /t 0 yields 

T c y/T 0 = 0.2489(14) , x 2 /d.o.f. = 1.5 . (8) 

The error bar here includes a rough estimate of systematic effects, encompassing the central 
values obtained by: (i) replacing “Wilson imp.” by “Wilson” or even the formerly excluded 
“clover” data; (ii) replacing the representation in eq. (4) through a) = Y2n=o a n(f3 — 

6.25) n ; (iii) carrying out the continuum extrapolation with a cubic fit; (iv) omitting f3 c 
corresponding to N T = 4 from the fit. The biggest deviations (T cy /t^ ~ 0.250) result either 
from using “clover” data which we assume to suffer from finite-volume effects, or from method 
(ii) which leads to x 2 larger by more than an order of magnitude in eq. (8). (An analysis 
based on data for t 0 /a 2 from previous literature can be found in ref. [8], is however subject 
to noticeably larger finite-volume effects than our current determination.) 

Comparing eq. (8) with eq. (6), we extract y/t^/r 0 = 0.3338(28), in perfect agreement with 
y/t^/r q = 0.3343(21) from refs. [15,35]. It is comforting to find a good agreement from a 
largely independent analysis. 


5. Conclusions 

In this note we have demonstrated that with modern resources and an opportune choice of 
an observable, the critical coupling j3 c of the Wilson plaquette action can be determined with 
<0.1% errors up to N T ~ 20 (cf. table 1). Subsequently, the critical temperature T c of pure 
SU(3) gauge theory could serve as a valid scale setting parameter for values of the Wilson 
coupling in the range 5.7 < f3 <6.8 (cf. table 1, from which the lattice spacing a is obtained 
as a = l/(N r T c ) if we simulate at the (3 C corresponding to N T ). Unfortunately these values 
are not large enough for scale setting on the very fine lattices (for instance N r = 48, (3 ~ 7.8) 
that are being used for studying transport observables close to the continuum limit [6]- [8]. 
Therefore “theoretical” quantities like r 0 and continue to be needed as intermediate 

3 For the sake of reproducibility of subsequent results we show more digits than are statistically significant. 



steps. On this point our study suggests that, with comparable numerical effort, employing 
y/t^ may yield more stable results than r 0 , however being assured that systematic errors are 
below the percent level remains a challenge for f3 > 6.4. If y/t^ is indeed used for scale setting, 
a conversion to T c can be carried out through eq. (8): y/t^T c = 0.2489(14). 

For various comparisons of lattice data with continuum perturbation theory, it would be 
very welcome to improve on our knowledge of y^A^g, whose uncertainty is currently an 
order of magnitude larger than that of y/t^T c . 4 Another issue worth further consideration is 
whether the method of sec. 2, which relied on the breaking of a discrete symmetry, could be 
generalized to the case of a continuous symmetry (such as a chiral symmetry). 
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